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Abstract: Nondestructive and noninvasive neutron assays are essential applications of 
neutron techniques. Neutron resonance transmission analysis (NRTA) is a powerful 
nondestructive method for investigating the elemental composition of an object. The 
back-streaming neutron line (Back-n) is a newly built time-of-flight facility at the China 
Spallation Neutron Source (CSNS) that provides neutrons in the eV to 300 MeV range. 
A feasibility study of the NRTA method for nuclide identification was conducted at the 
CSNS Back-n via two test experiments. The results demonstrate that it is feasible to 
identify different elements and isotopes in samples using the NRTA method at Back-n. 
This study reveals its potential future applications. 


Key words: CSNS, white neutron beam, NRTA, nuclide identification, nondestructive 
method. 


1. Introduction 


Neutrons are unique probes for analyzing the interior properties of materials in a 
nondestructive manner owing to their strong penetration power. Various neutron 
techniques based on the interactions between neutrons and materials have been 
successfully developed and applied in different fields. Neutron resonance analysis 
(NRA) is a technique that uses neutron resonance absorption to identify and quantify 
elements and isotopes in a sample [1—3]. Its basic principle is that the absorption of 


neutrons by nuclei as a function of neutron energy shows specific resonance peaks and 


dips based on which elemental composition of the material can be determined [4]. 
Neutron resonance transmission analysis (NRTA) is an NRA method based on the 
transmission of a neutron beam through a sample, in which the resonances are observed 
as dips as a function of neutron energy. 

The first demonstration of NRTA as a nondestructive method for determining the 
isotopic abundance of spent nuclear fuel pins was by Priesmeyer and Harz [3]. 7°°°8U 
and 7°°Pu contents were obtained in their study. Subsequently, NRTA for isotopic assay 
was used as a nondestructive technique to characterize special nuclear materials [5-7]. 
The results obtained using NRTA were consistent with those obtained using mass 
spectrometry and neutron activation methods [8]. NRTA has also been applied in other 
fields, such as the characterization of reference materials for nuclear data measurement 
[9] and cultural heritage objects and artifacts [10]. Moreover, NRTA has also been used 
with neutron imaging, known as neutron resonance transmission imaging, at the J- 
PARC facility [11], ISIS spallation neutron source [12], and KURRI-LINAC [13]. The 
NRTA method was investigated extensively by a group from EC-JRC Geel at the 
GELINA facility [14-16], which is a white neutron source that provides neutrons from 
10 to 20 MeV [17]. A white neutron beam covering the energy range from eV to MeV 
has significant advantages over thermal neutron beams because a broad energy range 
covers almost all the resonance regions of the elements, whereas the thermal beam is 
mainly used for heavy elements. The back-streaming neutron line (Back-n) at China 
Spallation Neutron Source (CSNS) is a newly built white neutron beam covering the 
range from 0.5 eV to 300 MeV [18,19]. The time-of-flight (TOF) technique was used 
at Back-n to accurately determine the neutron energy. The energy resolution at Back-n 
end station 2 (ES#2) was lower than ~2% for neutrons below 1 MeV. Because of these 
favorable beam properties, Back-n is suitable for the study and application of the NRTA 
technique. 

This paper describes the first test experiments using the NRTA method at the CSNS 
Back-n facility. The experimental results for two samples containing light or medium 
elements are presented. Data analysis shows that element/isotope identification using 


the NRTA method is feasible at Back-n. 


2. Experimental setup of NRTA at CSNS Back-n 


Located in Dongguan, Guangdong, the multi-disciplinary CSNS facility has been 
operational since 2018. It generates neutrons via a spallation reaction by impinging 1.6 
GeV protons onto a massive tungsten target. The nominal proton beam power of CSNS 
Phase I is 100 kW. In the forthcoming CSNS Phase II upgrading project, the beam 
power will be increased to 500 kW, and several more spectrometers [20,21] and 
experimental stations [22] will be built. Although most neutron beamlines at CSNS 
deliver thermal neutrons for neutron scattering applications, Back-n beamline delivers 
neutrons within a very broad energy range (from 0.5 eV to 300 MeV) mainly for nuclear 
data measurements [19,23,24] and nuclear technology applications [25,26]. 

Since its inception, the Back-n beam’s time characteristics have been a focal point 
of investigation [27] due to their crucial role in determining neutron energy and 
influencing energy resolution. The energy resolution of a TOF facility, considering the 


relativistic effect, is calculated as 
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where y is the Lorentz factor, v is the neutron velocity, and c is the speed of light. T and 
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L are the neutron flight time and flight path length, respectively, and AT and AL are 
their respective uncertainties. AT is mostly from the incident proton pulse width, whose 
full-width at half-maximum was 60 ns during the measurements in this study. AL is 
caused by the neutron moderation and scattering in the spallation target before they 
enter the neutron tube. The moderation process has been investigated via Monte Carlos 
simulations [28,29], and the uncertainties of the moderation distance (AL) at different 
energies were obtained. Table 1 lists the AL values at different energies and the energy 
resolutions calculated using Eq. (1), which indicates that the energy resolution is a 


function of the neutron energy. 


Table 1 Energy resolutions of CSNS Back-n at ES#2 
E, (eV) AT (ns) AL (cm) AE, /En 


1 60 12.2 3.14 x 10° 


10 60 13.6 3.50 x 10° 
10? 60 24.0 6.18 x 10° 
10° 60 20.2 5.24 x 10° 
104 60 18.0 5.10 x 10° 
10° 60 15.4 7.84 x 10° 
10° 60 8.8 2.15 x 10° 
10’ 60 10.0 6.82 x 10° 

~76 meters 

~55 meters 
End Station,1 (ES#1) End Station 2 (ES#2) 


Sample changer 


Neutrons Neutron flight tube 


Spallation 
target 


Shutter Collimator 1 Collimator 2 


Sample 
Fig. 1. Layout for neutron total cross-section measurements 

In this study, we implemented NRTA measurements using a neutron total cross- 
section spectrometer. This spectrometer comprised a multilayer fast fission chamber 
(FIXM) and a sample changer [30]. The FIXM is a neutron detector that uses 7*°U and 
°38LJ samples as neutron converters. Fig. 1 shows the typical setup for the neutron total 
cross-section measurements. The samples to be measured were placed on a sample 
changer at end station 1 (ES#1), and the detector was set up at ES#2. The beam spot at 
the sample position is typically turned to be ®60 mm, whereas that at the detector 
position is typically configured to be approximately ®40 mm. This configuration 
ensures that the neutron beam at ES#2 is entirely covered by a neutron converter with 
a diameter of D50 mm. Total cross-section measurements typically include sample-out 
and sample-in measurements. Sample-out measurements involve characterizing a 
neutron beam without any sample in the beam, whereas sample-in measurements 


involve placing the sample in the beam at ES#1. The long distance (approximately 20 


m) and the collimator between the sample and detector minimally suppress the 
multiscattering neutron background. The transmission spectrum is the ratio of the 
normalized neutron counts of the sample-in to the sample-out, from which the total 
cross-section can be determined. The neutron total cross-sections of ™C [30], “Li [31], 
Be, 77Al [32], Fe, ™'Pb, "Cr, and *°°Bi [33] have been measured using this setup at 
Back-n. 

The NRTA experiment at Back-n utilized the same setup as that used for the total 
cross-section measurements. Photographs of the sample changer and neutron detector 
are shown in Fig. 2. The resonances were observed as absorption dips in the 
transmission spectrum as a function of the neutron energy, from which the nuclides in 
the sample could be identified because the resonances of each nuclide are as unique as 
fingerprints. The energy positions of the resonance dips in the transmission spectrum 
provide qualitative information about the composition, whereas the profiles of the 
resonance dips present quantitative details about the composition. In our first NRTA 
experiments, we resolved the qualitative information about the samples. Quantification 


of the elemental composition is more complicated and will be our subsequent task. 


Fig. 2 Sample changer at Back-n ES#1 (left) and detector at ES#2 (right) 


3 Results and discussion 


The NRTA method is centered on measuring and analyzing neutron transmission 
spectra. The TOF technique was used to determine the neutron energy. Further details 
of the TOF method at Back-n can be found in our previous publications [18, 34]. 


Another crucial aspect of neutron energy determination is the double-bunch effect, in 


which two identical proton bunches have well-defined intervals (410 ns) in an incident 
pulse. Double-bunch unfolding is necessary for energy regions higher than 10 keV, as 
the 410 ns uncertainty is no longer negligible. An unfolding code based on Bayes’ 
theorem was developed to unfold the TOF spectrum [35,36]. The systematic uncertainty 
introduced by the unfolding process can be estimated using the unfolding code, which 
is highly dependent on statistical uncertainties. The uncertainty of the unfolding is 
typically 1.5—2.5 times that of the statistical uncertainty. 


3.1 Experiment using imitated lunar soil 


As shown in Fig. 3, an artificially simulated lunar soil sample was measured to 
verify the feasibility of the NRTA at Back-n. The sample was a cylinder with a diameter 
of 27 mm and a height of 21 mm, and comprised uniformly distributed elements. The 


details of its composition are listed in Table 2. 


Fig. 3 Photograph of imitated lunar soil sample 


Table 2 Composition of imitated lunar soil sample 
Composition SiO2 _AbO; CaO __Fe203 K20 MgO NaO TiO2 


Mass fraction (%) 49.6 13.8 9.6 14.8 0.48 7.85 2.07 0.922 

Fig. 4 compares the measured fission rate spectra of the sample-in and the sample- 
out. The sample-in spectrum (the red curve in Fig. 4) was attenuated by the scattering 
and absorption by the sample. Fig. 5 shows the transmission of the neutron beam 
through the sample, which is the ratio of sample-in to sample-out measurements. 
Resonance absorption dips were observed in the transmission spectrum and could be 


used to determine the existing isotopes in the sample. 
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Fig. 4 Comparison of sample-out and sample-in measurements of imitated lunar soil 
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Fig. 5 Transmission spectrum of imitated lunar soil 
Fig. 6 shows the transmission spectra in different energy regions to provide clear 
profiles of the absorption dips. All results are presented in 100 bins/decade (bpd), and 
the bin center is adopted as the energy of each bin. The resonances of sodium, aluminum, 
magnesium, iron, silicon, and oxygen were recognized. Table 3 presents a comparison 
of the measured resonance energy values and the corresponding resonance energy 
values in the ENDF/B-VIII.0 library [37]. The differences between the values of the 
measurements and library, and the energy resolution at each energy level are also 
presented. The differences at most points are comparable to the energy resolution, 


except at 88.11 keV, which is assumed to be the overlap of the resonances of *Mg and 


Fe at 83.45 and 84.84 keV, respectively. The overlap of the two resonances may 


increase the deviation between the measured and evaluated data. The differences vary 
with energy because, as discussed in Section 2, the energy resolution is a function of 
the neutron energy. Binning is another reason for the differences in the resonance 
energy positions between the measured and ENDF/B-VIII.0 library values. The energy 
bin center depends on the bin width, and smaller energy binning should provide a more 
accurate resonance energy; however, it increases the statistical uncertainties. This issue 
can be addressed using long beam-time measurements in the future. In this study, a 100 
bpd was used as the final option. In addition, although the unfolding method was used 
in a good order for the region higher than 10 keV, the double-bunch commission 


inevitably introduced systematic uncertainties. 
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Fig. 6 Expanded transmission spectra of imitated lunar soil in different energy regions 


Table 3 Resonance energy of imitated soil sample 


Resonance energy (keV) Energy resolution 
Isotope Difference (%) 

Measurement ENDF B-VIII.0 (%) 
Na 2.85 2.81 1.4 0.63 
Na 1.60 x 10° 1.60 x 10° 0 2.4 
27 Al 35.08 35.06 0.06 0.59 
Mg 88.11 83.45 5.58 0.75 
Mg 254.11 257.19 1.20 1.02 


Mg 431.55 433.08 0.35 1.29 


“Fe 7.67 7.82 1.96 0.80 


aoe 27.86 27.92 0.22 0.56 
Re 88.11 84.84 3.85 0.75 
?8Si 55.60 55.80 0.36 0.65 
160 431.55 434.00 0.57 1.29 
160 1.01 x 10° 1.00 x 10° 1.0 2.19 
160 1.30 x 10° 1.31 10° 0.77 2.31 
t 1.93 x 10° 1.90 x 10° 1.58 2.63 


Through qualitative analysis, most elements listed in Table 2 were identified, 
except for calcium, potassium, and titanium. Their absorption was not evident, probably 
because of the significant statistical fluctuations and their low abundance. In the above 
analysis, we only selected apparent dips and did not attempt to determine all the 


resonances of each nuclide. This is feasible for a qualitative study of the sample. 
3.2 Experiment using relic obtained from South China Sea 


The NRTA at Back-n also investigated an underwater relic found in the South 
China Sea. Fig. 7 shows a block of silicate and calcium carbonate embedded in many 
copper coins. It had an irregular shape, with an approximate size of 25 cm x 10 cm x 
10 cm. The experimental setup was the same as that used to measure the imitated lunar 
soil. However, this object contained nonuniformly distributed elements and had an 
irregular shape. An associated experiment using a CMOS camera was also conducted 
to capture transmission images of the sample using the total transmission and 


resonance-selected methods [38]. 


Fig. 7 Photographs of underwater relic found in South China Sea 


Fig. 8 displays the spectra for both sample-in and sample-out measurements. The 
sample-in measurements significantly attenuated because of the large sample size. Fig. 
9 shows the transmission of the neutron beam through the sample, from which the 
isotopes in the sample is examined. Fig. 10 shows zoomed-in graphs of the transmission 
spectra in different energy regions. The resonances of copper, calcium, sodium, silicon, 
and iron were observed. These recognized isotopes were consistent with our 
expectations. Cu and Cu originated from the coins in the object. *8Si, Ca, *°Fe, and 
160 originated from silicate and calcium carbonate, which are the main components of 


the object. The presence of Na is normal because sodium is abundant in the sea. 
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Fig. 8 Comparison of sample-out and sample-in measurement results of underwater relic 
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Fig. 9 Transmission spectrum of underwater relic 
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Fig. 10 Expanded transmission spectra of underwater relic in different energy regions 


We compared the measured resonance energy positions with corresponding values 
in the ENDF/B-VIII.0 library (Table 4). The energy binning was 100 bpd. The 
differences in the resonance energy positions between the measurement and the library 
were comparable to the energy resolution. In addition to the reasons mentioned in 
Section 3.1 for such differences, the irregular shape and massive volume of the sample 
could also increase the deviation between the measured and library values. 


Table 4 Resonance energy of underwater relic sample 


Resonance energy (keV) Energy resolution 
Isotope Difference (%) 
Measurement ENDF B-VIII.0 (%) 
"gi 55.60 55.80 0.36 0.65 
Na 2.85 2.81 1.4 0.63 
Fe 27.86 27.92 0.21 0.56 


56Fe 86.11 84.84 1.5 0.74 


Scu 0.582 0.578 0.7 0.59 


“Cu 2.07 2.05 0.98 0.61 
Cu 5.43 5.38 0.93 0.72 
“Cu 8.04 7.95 1.13 0.80 
Scu 2.54 2.54 0 0.62 
"Cu 3.94 3.92 0.5 0.67 
Scu 4.52 4.49 0.67 0.69 
&5Cu 6.53 6.45 12 0.75 
40a 133.36 133.02 0.23 0.83 
o 431.55 434.00 0.57 1.29 
o 1.01 x 10° 1.00 x 10° 1.0 2.19 
160 1.33 x 10° 1.31 x 10° 1.5 2.32 


4 Conclusions and prospect 


In 2018, CSNS introduced a new white-neutron beamline, Back-n. Its beam 
characteristics (high neutron flux, broad energy range, and good energy resolution) 
motivated us to launch a feasibility study of NRTA for nuclide identification. 
Experiments on two different samples, an imitated lunar soil sample and an underwater 
relic, were performed for the first time at Back-n. The existence of !°O, ?8Si, 77A1, *4>°Fe, 
24Mg, and Na in the simulated lunar soil sample and °-Cu, '!°O, 78Si, *°Fe, Ca, and 
°3Na in the underwater relic were confirmed using the NRTA method. A simple 
qualitative analysis demonstrated that CSNS Back-n could identify nuclides using the 
NRTA method. The beam characterization and results of this study demonstrates that 
the CSNS Back-n is a promising platform for nuclide identification. 

FIXM serves as the primary neutron detector for transmission measurements at 
Back-n. This helps in monitoring the flux based on 7°°U(n, f) and *8U(a, f) reactions. 
However, the resonances of 7*°U(n, f) cross-section below the ~keV region are 


somewhat intense, which might induce fluctuations in the transmission spectrum. 


Lithium glass can be a valuable compensation for the fission chamber in the low-energy 
region because of the large and smooth cross-section of the °Li(n, t) reaction [39]. A 
lithium-glass monitor is planned for future transmission measurements at Back-n, 
which is characterized by its quick response, high efficiency, and smooth transmission 
spectrum. However, because the lithium glass scintillator is sensitive to y-rays, n/y 
discrimination [40] and y-rays background must be investigated thoroughly. 

In this paper, we present a simple qualitative analysis focused on confirming the 
presence of specific isotopes. Accurate quantification analysis for resolving the 
elemental abundance in a sample is significantly more complicated. Quantification 
analysis is feasible by performing sophisticated resonance analysis, which must reflect 
the experimental effects, resolution function, sample properties, and detector 


characteristics. This is the subject of our subsequent study. 
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